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Abstract 

A simple many-fermion system in which there exists N identical fermions in a 
single spherical orbit with pairing interaction is treated by means of the time-dependent 
variational approach with a quasi-spin squeezed state with the aim of going beyond the 
time-dependent Hartree-Fock or Bogoliubov theory. It is shown that the ground state 
energy is reproduced well analytically in this approach. 
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§1. Introduction and Preliminary 



The time-dependent variational approach with the sii(2)-coherent state leads to the time- 
dependent Hartree-Fock (TDHF) or Bogoliubov (TDHB) theory for simple many-fermion 
systems governed by the dynamical sii(2)-algebra. In the previous paper, with the aim of 
going beyond the TDHF approximation, we have investigated the dynamics of the Lipkin 
model by means of the time-dependent variational approach with a quasi-spin squeezed state. 
The quasi-spin squeezed state is a possible extension of the sii(2)-coherent state. ^^'^^ Thus, 
the use of this extended trial state leads to an extended TDHF theory.^) In Ref.l), the 
role of the quantum effects included in the quasi-spin squeezed state was analyzed in the 
Lipkin model. Also, in Ref.5), the fermionic squeezed state was introduced for the 0(4) 
model with a schematic pairing plus quadrapole interaction and the role of quantum effects 
were analyzed. The use of the squeezed state in the time-dependent variational approach 
is one of possible candidates to go beyond the TDHF or TDHB approximation. However, 
in the previous works, we have only investigated the role of the fluctuations included in the 
squeezed state based on numerical results such as the ground state energy and dynamical 
motion. Thus, it is desirable to give an investigation based on an analytical results, especially, 
for the approximate reproduction of the ground state energy in dynamical viewpoint in order 
to clarify the role of quantum effects included in the quasi-spin squeezed state all the more. 

In this paper, paying an attention to the role of the quantum fluctuation included in the 
quasi-spin squeezed state and interesting to the dynamics of a system, we treat a simple 
many-fermion system in which there exists N identical fermions in a single spherical orbit 
with pairing interaction. The exact energy eigenvalue of the pairing model is known analyt- 
ically. Thus, we can compare the result obtained in our squeezed state approach with the 
exact one. The main aim of this paper is to give an analytical understanding for the role of 
the quantum effects included in the quasi-spin squeezed state. Especially, the dynamics of 
the variable representing quantum fluctuations is taken into account. As a result, it is shown 
that, under a certain condition, the ground state energy is reproduced well analytically, like 
expansion method. 

The single particle state is specified by a set of quantum number {j, m) , where j and m 
represent the magnitude of angular momentum of the single particle state and its projection 
to the 2;-axis, respectively. Let us start with the following Hamiltonian: 

mm m' 

where e and G represent the single particle energy and the force strength, respectively. The 
operators and c|„ are the fermion annihilation and creation operators with the quan- 
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turn number m, which obey the anti-commutation relations: { Cm , cl^/ } = ^mm' and 
{ Cm 1 Cm' } = { 1 c^m' } ~ ^- introduce the following operators: 

'^+ = ^E^-4> = ^o = ^E4cm-^?) , (1-2) 

m m m 

where Cfh — {—y~'^C-m and O represents the half of the degeneracy: f2 = j + 1/2. These 
operators compose the sw (2) -algebra: 

[^+,^_]=2^o, [S-c^il^i^i . (1-3) 

Thus, these operators are called the quasi-spin operators.^-*' Then, the Hamiltonian (IT) 
can be rewritten in terms of the quasi-spin operators as 

H = 2e(^o + Sj) - GS+S- ^ eN - GS+S- , (T4) 

S, = n/2 (= (j + l/2)/2) , 

where N represents the number operator: 

N=J2^lnCm. (1-5) 

m 

As is well known, the eigenstatcs and eigenvalues for this Hamiltonian are easily obtained. 
Then, the ground state energy can be obtained by setting the seniority number being zero:*^ 

£„ = .^-lG/V«(2-§ + i) . (1.6) 

Next, we review the coherent state approach to this pairing model, which is identical 
with the BCS approximation to the pairing model consisting of the single energy level. The 
sw(2)-coherent state is given as 

10(a)) =exp(/^+-r^_)|0) , 

|0) = |5 = Sj, So = -Sj) , (^_|0) = 0) . (1-7) 
We impose the canonicity condition :^^ 

(0(a)|^|0(a)) = IC , (0(«)l^|0(«)) = • (1-8) 
A possible solution of the above canonicity condition is presented as 

e = v^^sin|/|, r^v^j^sinl/l . (1-9) 
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Thus, the expectation values of Hamiltonian H and number operator are easily calculated 
and given as 

{^{a)\H\^{a)) = 6 . 2|^|^ - G (28^? ' {} ' ^) '^'') ^ ^ ' 
{^{a)\N\4^{a)) ^ 2\i\' = N . (1-10) 

If total particle number conserves, that is, N ^constant, then, the energy expectation value 
E is obtained as a function of N : 

1 / N N \ 

E = eN--GNn{2-- + ^^ . (Ml) 

Compared the approximate ground state energy in (1-11) with the exact energy eigenvalue 
in (1-6), the last terms in the parenthesis of the right-hand side are not identical, which has 
order of 1/^2 if the order of the magnitudes of N and f2 is comparable. The main aim of this 
paper is to show that the last term, 2/J7 in (1-6), can be recovered by taking into account 
of the dynamics in the time-dependent variational approach with the quasi-spin squeezed 
state. 

§2. Quasi-spin squeezed state for pairing model 

In this section, the quasi-spin squeezed state is introduced following to Ref.2). First, we 
introduce the following operators: 

it = , A = , TV = 2{So + Sj) , (2-1) 

where N is identical with the number operator (1-5). Then, the commutation relations can 
be expressed as 

[A, it] ^ 1 _ ^ , [TV, i] = -2i , [TV, it] = 2it . (2-2) 

Using the boson-like operator i^^ sM(2)-coherent state in (1-7) can be recast into 

10(a)) = -^=l=exp(ait)|0) , ^a*a) = (^1 + , (2-3) 

where a is related to / in (1-7) as a = ^/2Sj ■ (//|/|) ■ tan |/|. The state |0(q;)) in (2-3) is a 
vacuum state for the Bogoliubov-transformed operator : 

am = Ucm - V{-y-^cl^ , a^|0(a)) = . (2-4) 



4 



The coefficients U and V are given as 

, ^ " ^ ^ c/2 + |y|2 = i. (2-5) 

V'l + a*a/(25'j) ^^25^ ^1 + a*a/(25'j) 

Of course, and al^ are fermion annihilation and creation operators and the anti-commutation 
relations are satisfied. By using the above Bogoliubov-transformed operators, we introduce 
the following operators: 



V J m V J m 

Then, the state satisfies 



= M\<t>{a)) = . (2-7) 
Further, the commutation relations are as follows: 

[B, St] = 1 - ^ , [M, B] = -2S , [M, St] = 2St . (2-8) 

The quasi-spin squeezed state can be constructed on the sii(2)-coherent state |0(q;)) by using 
the above boson-like operator W as is similar to the ordinary boson squeezed state: 



*(/'-/')-i+|^n(i-4)w)'. (2-9) 

We call the state \'4>{oi,l3)) the quasi-spin squeezed state. 

We can easily calculate the expectation values for various operators with respect to the 
quasi-spin squeezed state. The expectation values can be expressed in terms of the canonical 
variables which are introduced through the canonicity conditions. The same results derived 
in the following are originally given in the Lipkin model in Ref.lO) at the first time and 
these results were used in Ref.l) to analyze the effects of quantum fiuctuations in the Lipkin 
model. For the quasi-spin squeezed state \'^{a, (3)) in (2-9), the following expression is useful 

(2.10) 
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where dz = d/dz. The canonicity conditions are imposed in order to introduce the sets of 
canonical variables {X,X*) and {Y,Y*) as follows : 

{^/;{a,(3)\dYma,(3)) = , {^|;{a, P)\dy.ma, P)) = -^F . (2-11) 

Possible solutions for X and Y are obtained as 

X , X* 

a — — = , a — 



I _ _ 4y*y /-^ _ X*X _ 4,Y*Y 

Y Y* 
/3 = , , /3* = , , (2-12) 

^/K{Y*Y) VK{Y*Y) ^ ^ 

where K{Y*Y) is introduced and satisfies the relation 

K{Y*Y)^{Y*Y/K) = ^\Y*Y/K) . (2-13) 

The expectation values for B, B\ M and the products of these operators are easily obtained 
and are expressed in terms of the canonical variables as follows : 

{^{a,(3)\B\i,{a,f3)) = {i^ia, f3)\B^\^P{a, f3)) = , 

{ipia, /3)|M|^(«, (3)) = 4Y*Y , (2-14a) 
{ij{a,P)\B'\^{a,P)) = 2Y^K{Y*Y) , 



(V'(a,/?)|StXa,/?)) = 2Y*^K{Y*Y) , (2-14b) 
{^l;(a,(3)\B^B\^j;(a,(3)) = 2(1 - l/(2Sj))Y*Y - ^{Y*Yf ■ L{Y*Y) , 

{ij{a,P)\BB^^{a,P)) = {^{a, P)\B^ B\^{a, P)) + (1 - 2Y*Y/Sj) , (2-14c) 

{ip{a, (3)\M^\ip{a, (3)) = 16r*r (1 + Y*Y ■ L{Y*Y)) , (2-14d) 

where L{Y*Y) is defined and satisfies 



K{Y*Yf . L{Y*Y) = . (2-14e) 

By using the relations between the original variables a and P and the canonical variables 
X and y, the coefficients of the Bogohubov transformation (2-5), U and V, are expressed as 



^ _ X'X _ 4Y*Y 

U = ^ , . V= ,„„ , . (2-15) 



y 
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Then, the operators A, and N, which are related to the quasi-spin operators S-, and 
5*0, respectively, in (2-1), can be expressed as 

A = I 1 - ^1 - V'^B'^ + U'^B , 



it ^ ,j2SjUV* 1 



M 
2Sj 



+ U^B^ - V*^B , 



M 



N = ASjV*V I 1 - ^ 1 + V^jU{VB^ + V*B) + M . 



(2-16) 



Thus, the expectation values for A, A^, N and the products of these operators are easily 
obtained and are expressed in terms of the canonical variables. For example, from (2-1), the 
expectation values of quasi-spin operators are derived as 



{ilj{a, P)\S+\ijj{a, (3)) = X*^j2Sj - X*X - 4Y*Y , 
{'ijj{a,P)\S_\'ilj{a,f3)) = ■s/2Sj - X*X - AY*Y X , 
{ij{a,P)\So\ip{a,P)) = X*X + 2Y*Y - Sj 

and also the expectation value of the number operator in (1-5) is calculated as 



(2-17) 



N = {iPia, (3)\N\ip{a, (5)) = 2X*X + AY*Y . 



(2-18) 



The above expressions in (2-17) correspond to the Holstein-Primakoff boson realization for 
the sM(2)-algebra such as 5*+ = X* ^2Sj — X*X and so on. Thus, we can conclude that the 
variable \Y\'^ represents the quantum effect. 

The model Hamiltonian (1-4) can be expressed in terms of the fermion number operator 
N and the boson-like operators A and A"^ as 



H^eN - 2SjGA'^A 



(2-19) 



Thus, the expectation value of this Hamiltonian is easily obtained. We denote it as ifgn : 



sq 



(V^K/3)|^|V^(«,/?)) 



= e{2nx + 4ny) - 2GSj |a' + 



- — ^/^Xycos{2ex -9y) + 2 



n 



X 



4ny) 

xy 



- -^Xy [ny + nlL - nl] 



S] 



1 - 



2S, 



ny — 



riy ■ L 



Si 



(2-20) 
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where we introduce the action-angle variables instead of {X,X*) and {Y,Y*) as 

X = V^e-''^ , X* = v^e^^- , 

Y = v^e-^^^ , y* = Vn^e'^'' (2-21) 
and A" and y are defined as 



1 - 



The dynamics of this system can be investigated approximately by determining the time- 
dependence of the canonical variables {X,X*) and {Y,Y*) or {nx,Ox) and (ny,^y). The 
time-dependence of these canonical variables is derived from the time-dependent variational 
principle : 

S J{tp{a, P) \idt - H\ip{a, P))dt = . (2-23) 
§3. Time evolution of veiriational state 

In the sM(2)-coherent state approximation, the expectation value of the Hamiltonian is 
calculated as 

//eh = = 2enx - 2GSj jn^ - || + ^| • (3-1) 

In this approximation, namely, usual time-dependent Hartrce-Bogoliubov approximation, 
the canonical equations of motion derived from the time-dependent variational principle 
have the following forms: 



dnx """'V Sj 2S] 

The solutions of the above equations of motion are easily obtained as 



Ox it) 



xo 



nxit) — no (constant) . (3-3) 

On the other hand, in the quasi-spin squeezed state approximation, the equations of 
motion derived from (2-23) are written as 

^ a — 
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S,{2Sj-AnY)^ S] dnx 



■ 1^ cos(2^x - ) - 2 ■ 1^ ^ ■ L 



dn 



X 



dnx Sj 



nx 



OH. 



sq 



de 



X 



(34) 



sq 



dny 

4e - 2GSj 



dX , n\-y 



+ 



4 



2S^ ^ S] 



dn 



dny 
2ny] 



■ y ■ [ny + UyL — ny] 



T o 7V' 

-——={K + nyj—)Xy cos{2dx - Oy) 
Sjy/nyK any 



— — \pnyK 



dn 



Y 



dn 



Y 



cos(2^x - Oy) 



dny dny 

xy 



1 - 



1 



1 



1 

25^ 
1 

s. 



ny 
ny - -L 

or 2 dL 

2nyL + ny- 



dH. 



ny 



sq 



de 



-AG^/nyKXy sm{2ex - 6 



3 / "^0 V. 

y) 



dny 



y 



(3-5) 





s 



time 

Fig. 1. The time evolution of nx(i) and ny(t) 
is plotted in the case of the quasi-spin 
squeezed state approach. The parame- 
ters are taken as G = 1.2, e = 1.0 and 
17 = 2N = 8. The initial values are 
nx{t = 0) = 1.3, 0x{t = 0) = 0.0, 
ny{t = 0) = 0.35 and 6'y(t = 0) = 0.0. 



-35 I : : : : : : : : : 1 

2e+0OB 4e+00B 6e+006 ee+006 1e+007 1.2e+007 l.-le^OO? 1.6e+007 1.86+007 2e+007 

time 

Fig. 2. The time evolution of 0x{t) and Oy(t) 
is plotted in the case of the quasi-spin 
squeezed state approach. The parame- 
ters are taken as G = 1.2, e = 1.0 and 
17 = 2N = 8. The initial values are 
nx{t = 0) = 1.3, dx{t = 0) = 0.0, 
ny{t = 0) = 0.35 and Oyit = 0) = 0.0. 
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Here, the number conservation is satisfied: 

TV = 2nx + 4ny = . (3-6) 

The time evolution of nx{t) and ny(t) in Fig.l and Ox{t) and in Fig.2 is plotted 

with appropriate initial conditions. The parameters used here are G — 1.2, e = 1.0 and 
Q = 2N = 8. In the SM(2)-coherent state approach, nx is constant of motion. However, 
quasi-spin squeezed state approach, nx and ny oscillate with antiphase, while the total 
particle number is conserved. 

§4. Dynamical approach to the ground state energy 

Hereafter, we assume that |y p (= ny) <^ 1 because |y p means the quantum fluctuations. 
Then, K{Y*Y) and L{Y*Y) defined in (2-13) and (2-14e), respectively, can be evaluated by 
the expansion with respect to As a result, we obtain 

1 / 2 \ / S 

L{Y*Y) = r3 1- 



1 V '^^jJ\ '^^j 

' ^' ^)(^-^)(^-^]y'y+-- (4-1) 



(1 — 2^)^ V SSj/V 25*^ 

Then, the expectation values for the Hamiltonian, the time-derivative and the number op- 
erator can be expressed as 



^sq = e{2nx + ^ny) - G 



2S^nx-n\ + ^^ 



-2V2y'l-^ - nxV^cos(2^x - Oy) 

(10 2 8 \ 

2S^ - 1 - 4nx + ^nx + - j^^xj ny + ©(nj/^) 

(4-2a) 

{iP{a,(5)\%dt\^{a,(3)) = {nxdx+nyOy) , (4-2b) 
N = (^(a,/3)|7V|^(a,/3)) = 2nx + 4ny . (4-2c) 

Prom the time-dependent variational principle (2-23) or (3-4) and (3-5), the following equa- 
tions of motion are derived under the above-mentioned approximation : 



ex^^^2e-G-2 
dnx 



■V2^1-^(l - ^)v^cos(2^x - Oy) 
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+4 



sq 



de 



X 



(4-3a) 
(4-3b) 



sq 



dn 



Ae-G 



Y 



2Sj 



-V2jl-J^(l-^]^ cos(2^x - 9y) 



2S, 



+2 



1 



-1 + 2,5,- - 4nx + TT + TT^x - 7;^"' 



s. 



2"'X 



sq 



G ■ 2v^^l-^(^l - ||-^nxy^sm(2^x - By) . 



(4-3c) 
(4.3d) 



It is found from (4-3b) and (4-3d) that the total fermion number N in (4-2c) is also conserved 
in this approximation, that is, 

TV = 2nx + 4nr = . (4-4) 

It should be noted here that nx < N/2 from (4-2c) and N < Qj2 — Sj. Thus, the 
inequality 1 — nx/2Sj > is obtained. From the approximated energy expectation value 
(4-2a) for G > 0, the energy minimum is then obtained in the case cos(2^x — ^y) = — 1, 
namely, 

9y^29x + n. (4-5) 

Since the energy is minimal in the ground state, the relation (4-5) should be satisfied at any 
time. In order to assure the above-mentioned situation, the following consistency condition 
should be obeyed : 

Oy = 2ex . (4-6) 

Thus, from the equations of motion (4-3a) and (4-3c), under the approximation of small 
ny, the consistency condition (4-6) gives the following expression of ny in the lowest order 
approximation of ny : 



1-^ 



2Si 



2\/2(l-i) 



1 + 



1 



2nx{l-m{l-4-) 



-1 



(4-7) 



Thus, by substituting (4-5) and -y/ny in (4-7) under the lowest order approximation of ny 
into the energy expectation value (4-2a), and by performing the approximation of large N 
or large f2{— 2Sj) approximation, we obtain the ground state energy as 

^eN- ^GNf2 (^2-^ + ^ + 0{l/{Nf2), l/f2', l/N')^ . (4-8) 
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This result reproduces the exact energy eigenvalue (1-6) by neglecting the higher order term 
of 1/{NQ)^ l/Q"^ and for large and i? limit. Thus, the quasi-spin squeezed state 

presents a good approximation in the time-dependent variational approach to the pairing 
model. In this approach, the existence of the rotational motion in the phase space consisting 
of {nx,dx',nY,OY) plays the important role. The angle variables for rotational motion, 6x 
and 6y, are consistently changed in (4-6). This consistency condition is essential to reproduce 
the exact energy for the ground state under the large N and O limit. The approximation 
corresponds to so-called large approximation. In general, it is known that the large 
expansion at zero temperature corresponds to h expansion. In this sense, the time-dependent 
variational approach with the quasi-spin squeezed state gives the approximation including 
the higher order quantum fluctuations than h if any expansion is not applied. 

§5. Summary 

In this paper, it has been shown that the exact ground state energy for the pairing model 
can be well recovered by using the time-dependent variational approach with the quasi-spin 
squeezed state. For this purpose, we treated the sii(2)-algebraic model because its eigenvalue 
is known analytically. As a result, by taking into account of the dynamics in our quasi-spin 
squeezed state approach, the exact ground state energy can be reproduced up to the order 
of l/i7 under the small \Y\'^ expansion. Of course, the time evolution of a system governed 
by the pairing model Hamiltonian can be also investigated as is similar to that of the Lipkin 
model developed in Ref.l). However, we do not repeat it because the result is almost same 
as that of the Lipkin model, which was reported in Ref.l). 
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